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Abstract 

The kinetic evolution of the Orszag-Tang vortex is studied using coUisionless hybrid simula- 
tions. In magnetohydrodynamics this configuration leads rapidly to broadband turbulence. At 
small scales, differences from magnetohydrodynamics arise, as energy dissipates into heat almost 
exclusively through the magnetic field. A key result is that protons are heated preferentially in 
the plane perpendicular to the mean magnetic field, creating a proton temperature anisotropy as 
is observed in the corona and solar wind. 
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The dissipation of turbulent energy in plasmas plays a critical role in understanding 
coronal heating and the acceleration of the solar wind 111, turbulence in the interplanetary 

n n 

medium [2[, energy storage and release in the magnetosphere [3|, and in a variety of other 
plasma and astrophysical contexts [4]. By "dissipation," we mean the conversion of fluid 
scale energy irreversibly into kinetic degrees of freedom. A key observational clue about 
the nature of this dissipation is the substantial heating of protons, often preferentially in 
the plane perpendicular to the mean magnetic field[l|,l3], both in the solar wind (plasma 
f3 = (thermal speed/ Alfven speed) ^ ~ 1) and in the corona {(3 << 1). This heating has a 
variety of potential sources, including shocks and wave-particle interactions involving non 
thermal distributions such as pickup ions, but the ubiquity of broadband Kolmogoroff- 
like fluctuations suggests that kinetic absorption of fluid energy at or beyond the high 
wavenumber end of the inert ial range plays an important role. 

A fundamental demonstration of turbulent anisotropic heating is needed as a first step 
towards a basic physics understanding of the dissipation processes that heat the solar corona 
and solar wind. Magnetohydrodynamics (MHD) is a very useful plasma model which gener- 
ally employs small but nonzero viscosity and/or resistivity, but this is not easily justified for 
collisionless systems where the mean free path is comparable to the system size (as in most 
of the corona and solar wind). More sophisticated attempts to numerically model plasma 
dissipation, e.g., by employing hyperresistivity, hyperviscosity or indirectly by including the 
Hall or finite Larmor radius effects (see [s]), include only selected approximations to kinetic 
effects, and therefore do not include the wider range of mechanisms available to the kinetic 
plasma. 

In this paper, we report a demonstration of turbulent anisotropic proton heating in the 
Orszag-Tang vortex [tI] using a hybrid simulation, which includes all proton kinetic effects. 
The hybrid simulation results are very similar at large length scales to MHD simulations of 
the same system, but show significant differences at small scales where kinetic effects are 
important. Analysis of the hybrid results show that energy is dissipated into proton heating 
almost exclusively through the magnetic field and not through the proton bulk velocity. 
The proton heating occurs preferentially in the plane perpendicular to the mean magnetic 
field. These simulations, to our knowledge, are the first self consistent demonstration of 
turbulent anisotropic proton heating. Finally, effective transport coefficients from the hy- 
brid simulations are calculated, showing that the approximation of constant resistivity t] is 
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potentially reasonable (although it cannot reproduce the proton temperature anisotropy), 
but a constant viscosity v is untenable. 

Q 

The Orszag-Tang vortex[7|] is a well studied MHD initial configuration given by 

B = — sin ?/ x+sin2a; y + -Bg z; v = — siny x+sina; y (1) 

with B the magnetic field [Bg a uniform guide field) and v the proton bulk velocity in 
normalized units described later. This configuration leads immediately to strong nonlinear 
couplings, producing cascade-like activity that might reasonably approximate the highest 
wavenumber decade of the inertial range. These couplings, which are dominantly local in 
wavenumber, in turn drive the dissipation range. The physics of the Orszag-Tang vortex 
has been previously studied using incompressible [3] and compressible [8;] MHD simulations. 
Its robust production of nonlinear activity is a motivation for its frequent use in validating 



numerical schemes (see e.g. IqI). 

Simulating kinetic dissipation is difficult and computationally expensive due to the re- 
quirement of treating a wide range of length scales. By choosing a computational domain 
with approximately one decade of scale in the MHD range and another in the kinetic range, 
we can study the conversion of strongly driven MHD fiuctuations into kinetic motions. An 
antecedent of the present study compared global behavior of hybrid and Hall MHD simula- 
tions lOj, but included a mean in-plane magnetic field while not adequately resolving the 



proton inertial length. 



We use the hybrid code PSD in 2.5D, which models protons as individual particles 



and electrons as a fiuid and evolves 

— = Vx(vxB)-6„Vx (-xB'j (3) 

where J = V x B is the current density, en = c/{LQUJpi) is the normalized proton inertial 
length, rrie and rrij are the electron and proton masses, Xj and Vj are the positions and 
velocities of the individual protons, and v is the proton bulk fiow speed. Length is normalized 
to Lq, velocities to Vq = Bo/(ATTmnoY^'^ , time to to = Lq/Vq, and temperature to -Bg/(47rno). 
The average density is no, and Bq is the root mean square in-plane magnetic field. The 
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FIG. 1: (Color online) Current density with magnetic flux contours at t=1.96 in (a) fluid and (b) 
hybrid simulations. 

magnetic field B is determined from B' using the multigrid method. The code assumes 
quasi-neutrality. The electron temperature is zero and is not updated. 

Hybrid simulations are ideally suited for exploring dissipation and proton heating in 
coUisionless plasmas because they include a complete kinetic description of protons. Due to 
the finite temperature of the protons, kinetic Alfven waves are present in this set of hybrid 
simulations ({12!. and references therein), as well as parallel proton bulk flows. 

The simulation domain is a square box of side length 2 7r x 2 7r with 512 x 512 grid 
points. About 2.6 million protons are loaded with an initial Maxwellian distribution having 
a uniform temperature = 8, and en = 27r/25.6 and rrie = 0.04mj. No artificial dissipation 
is present other than grid scale dissipation. Choosing a guide field Bg = 5 (total /3 = 
2nT/B'^ ^ 0.62) reduces the system compressibility. Incompressibility is further promoted 
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FIG. 2: Hybrid and MHD comparison: (a) Magnetic energy Eb, fluid flow energy E^, their sum, 
the change in thermal energy AEth, and total energy Etot vs. time, (b) Flow enstrophy < > 
and magnetic enstrophy < > vs. time. 

by adding perturbations to the background density uq that enforce 9( V ■ v)/(9t = at t = 0. 
Simulations without the added perturbation show only small differences. 

The hybrid simulation results are compared to those of a compressible 2.5 D MHD version 
of the code F3D 3 

with constant and uniform resistivity rj = 0.0048, zero viscosity z/, and 
ratio of specific heats 7 = 5/3. (We motivate values for rj and u later.) In both cases, the 
magnetic islands initially centered on the midplane {y = tt) begin a clockwise rotation. The 
initial velocity profile shears the magnetic islands until t ~ 2 as the islands approach and 
undergo a brief period of magnetic reconnection from t ~ 2 — 4. After t ~ 4, the system 
is dominated by strong turbulence. A comparison of out-of-plane current density Jz and 
magnetic field lines at t = 1.96 is shown in Fig. [H The hybrid and MHD results show strong 
similarities at large scales, but significant differences at small scales where kinetic effects 
become important. 

We can quantify the differences by comparing energy and dissipation budgets, where flow 
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energy Ey = (p|v^|/2), magnetic energy Eb = (|Bp/2), thermal energy Eth (total proton 
kinetic energy minus flow energy), and total energy Etot- (• • •) denotes a volume average. 
Grid scale fluctuations in the hybrid data are smoothed using a standard local, weighted 
iterative averaging. 

Figure [21(a) shows E^jEb, their sum, AEfh (where A means the change since t = 0) 
and Etot, as a function of time from the hybrid and MHD simulations, with Eb and Etot 
shifted down by a constant for convenience. Note that Etot changes very little over the 
course of the hybrid run, demonstrating good numerical energy conservation. During the 
initial phase {t < 2), bulk flow energy is converted strongly into magnetic energy as fleld 
lines are stretched, but with little proton heating. The magnetic energy converts back to 
flow energy (with some heating) during the reconnection event (t ~ 2 — 4). Until t ~ 4, the 
energetics of the hybrid and MHD results are very similar. However, in the turbulent phase 
(t > 4), the hybrid and MHD codes show signiflcant differences, and more dissipation occurs 
in the MHD case. Notably, the proton thermal energy increases monotonically during the 
turbulent phase in the hybrid simulation, even without explicit dissipation. 

In MHD, mean square gradients of v and B are proportional to the energy dissipation 
rate. Although the hybrid code lacks explicit dissipation, it is instructive to compare in 
Fig. the out-of-plane enstrophy (mean square vorticity) {u"^) = (|z ■ (V x v)p) and out- 
of-plane magnetic enstrophy (mean square current density) (cj^) = (|z- (V x B)p) = (J^) in 
the hybrid and MHD simulations. At early time (t < 4), the hybrid and MHD enstrophies 
peak at about same time, but their magnitudes are very different, indicating that the length 
scales in the hybrid case are larger, probably due to flnite Larmor radius effects. During 
the turbulent phase (t > 4), the enstrophies continue to be different but the magnetic 
enstrophies are surprisingly similar. This suggests that the kinetic dissipation may have 
some resemblance to a classical resistivity, which we revisit later. Note that the value of 
enstrophy in the hybrid case is necessarily sensitive to the averaging that deflnes the fluid 
scales. 

To quantify the dissipation, consider the flow of magnetic energy in the system. Dotting 
the induction equation [Eq. Q] with B, averaging over space, and integrating over time 
gives 



where A refers to the change since t = 0. The v ■ ( J x B) term, called AE^b, is the exchange 
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FIG. 3: (a) AEb- change of Eb in the hybrid run; AE^b- exchange between Ejj and Eb', AEe'. 
electron kinetic energy; Vb- sum of these, total Eb dissipated, (b) and Vb are cumulative 
dissipation through bulk flow and magnetic channels, Dtot their sum, AEth change in thermal 
energy, (c) Parallel and perpendicular proton temperatures vs. time. 

in energy between bulk flow and the magnetic field and the d"^ term is essentially the electron 
kinetic energy. We define Vb &s the cumulative energy dissipated from the magnetic channel. 
Included in this term are non-MHD dissipative processes and grid scale dissipation. The 
first three terms are calculable from the simulations, so we may compute Vb- 

Fig. M^a) shows that for t < 4, Eb increases due to input from E^, then decreases as 
it is converted back. In the turbulent phase (t > 4), there continues to be conversion into 
magnetic energy. However, Eb decreases during this period, and both energy converted from 
the flow, and an approximately equal amount previously stored in the field, is absorbed in 
the magnetic dissipation Vb- 

A similar analysis can be done for energy flow in the bulk flow channel. Dotting the 
MHD momentum equation with v and integrating over time and space gives 




(6) 
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FIG. 4: Effective (a) resistivity r] and (b) viscosity vs. time. 

where P„ is the cumulative energy converted into heat from the flow channel through non- 
fluid effects and compression. Fig. [3](b) shows (dashed line), Vb (dot-dashed), Vtot = 
Vb + 'Dy (solid), and ^Eth (dotted). In the turbulent phase (t > 4), there is essentially no 
energy dissipated through the flow channel (P^, is constant). The energy dissipated through 
the magnetic fleld T>b traces very closely the increase in thermal energy of the protons AEth, 
so the main source of dissipation in this system is through magnetic interactions. The small 
departure between the two is accounted for by numerical heating, seen in the change in total 
energy Etot [see Fig. |2](a)]; this is only about 10 % of the total dissipated magnetic energy 
It is reasonable that little dissipation occurs through the flow channel. In the MHD regime 
(low wavenumber k), the dynamics of the hybrid simulation are at most weakly compressible, 
so the majority of energy is in oblique Alfven waves that are weakly damped [l^. As energy 
cascades to smaller scales, the proton gyroradius is reached before the dissipation scale. 
Below the proton gyroradius, the ions decouple from the magnetic fleld and only weakly 
participate with the non-MHD waves in this region. Consequently, the Alfven ratio E^/Eb 
goes to zero in the kinetic regime as evidenced by the structure of kinetic Alfven and whistler 
waves. 

A central result of this study is that the dissipated magnetic energy preferentially heats 
the protons perpendicular to the mean magnetic fleld, as shown in Fig. [3](c). The perpen- 
dicular and parallel temperatures T± and T\\ are calculated relative to the guide fleld. In 
the turbulent phase (t > 4), T± increases monotonically, while T\\ remains relatively steady. 
The relative anisotropy is small because the available magnetic free energy in the system 
(from and By) is small compared to the proton temperature, i.e. j3± 3> 1. If the heating 
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mechanism in this study generahzes to [3\_ ~ 1, large anisotropies will develop. The per- 
pendicular heating occurs without any obvious connection to classical cyclotron resonances, 
as the latter generally are construed 16[ to involve gyroresonance with waves propagating 
parallel to a background field (along the invariant direction in this study). Resonance can 
also involve wave frequencies near the cyclotron frequency; however, dominance in this sim- 
ulation of incompressible modes and relatively low frequency kinetic Alfven waves make the 
connection to standard cyclotron resonance 16| uncertain. 

By assuming the classical functional forms for the dissipation rates rj{Jl) and i^{uj'^), 
we can compute effective transport coefficients rjes and from the hybrid simulations. 
Figure m shows r^efr = {dV b / dt) j {Jl) and v^a = (dVy/dt)/ versus time. Surprisingly, 
the spatially averaged rj^s is fairly constant in time, as is often assumed in MHD models. 
The mean value is rjes = 0.0048, corresponding to a magnetic Reynolds number of Ses = 
AttVoL / rjcsc^ ~ 1308, which is the value used in the MHD simulation. In classical turbulence 
theory, the length scale at which dissipation occurs is related to the Reynolds number 
and energy containing scale L through Ses ~ {L/XdY^^ -^^^ hybrid simulation, 

\d ~ 0.029, which is of the order of the electron skin depth, c/uope ~ 0.049. The spatially 
averaged z/eg, on the other hand, shows oscillations much larger than the mean, calling into 
question the assumption of a non-zero viscosity assumed in many MHD models. 

The effective resistivity being fairly constant does not imply that the dissipation is of the 
form rjJ^. The MHD simulations performed with this ?7efr show more dissipation of B than 
the hybrid simulations, and cannot reproduce the preferential heating of T^. Future studies 
will investigate the spatial dependence of T^es, and its dependence on electron/proton inertial 
lengths, as well as system size. We will also examine the effect of reducing the guide field and 
including electron pressure in Ohm's law. The physical mechanism which converts magnetic 
energy into proton heat remains an open question and is under investigation. Potential 
explanations are wave-particle interactions such as Landau damping. 

This work is supported by NSF ATM-0539995, NASA NNG05GM98G and NNX-08AI47G 
(Heliophysics Theory program), and NERSC. The authors thank S. Servidio, S.P. Gary, and 
B. N. Rogers for enlightening conversations. 
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